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We demonstrate that particles confined to two dimensions (2d) and subjected to a one- 
dimensional (Id) periodic potential exhibit a rich phase diagram, with both "locked floating solids" 
and smectic phases. The resulting phases and phase transitions are studied as a function of tem- 
perature and potential strength. We flnd reentrant melting as a function of the potential strength. 
Our results lead to universal predictions consistent with recent experiments on 2d colloids in the 
presence of a laser-induced Id periodic potential. 
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Landau's order parameter expansion predicts that the 
direct transition from a sohd to a hquid should always be 
first order 0|. In two dimensions (2d), however, fluctua- 
tions can suppress the transition temperature so far be- 
low its mean-field value that order parameter amplitude 
fluctuations (except in the form of topological defects) 
play no role, and Landau's mean-field analysis is quali- 
tatively wrong. In this case a two-stage melting process 
mediated by the unbinding of dislocations ||,^ and discli- 
nations ||^ provides an alternative scenario with two suc- 
cessive continuous phase transitions with an intermediate 
hexatic phase instead of a single direct first order transi- 
tion. A periodic embedding medium for the 2d solid (e.g. 
a crystal substrate or a laser potential) leads to commen- 
surability effects and engenders an enormous variety of 
interesting phenomena 

Colloids confined between smooth glass plates pro- 
vide an ideal model system for experimental studies of 
2d melting. In this system individual particles can be 
imaged, allowing a direct observation of topological de- 
fects and measurement of real-space correlation func- 
tions. Murray etal. |^] and Zahn etal. Q have given 
strong experimental evidence for a two-stage melting 
mechanism in such systems (because dislocations relax 
slowly, longequilibration times are required |^). Chowd- 
hury etal. jj] studied the behavior of a 2d colloidal sus- 
pension subjected to a Id periodic potential provided by 
the standing wave pattern of two interfering laser beams. 
They observed a phenomenon called "light induced freez- 
ing" , a solidification driven by a commensurate laser po- 
tential. Our research is motivated by the recent work of 
Wei etal. who reinvestigated this phenomenon and 
discovered reentrant melting upon gradually increasing 
the amplitude of the laser potential even further. 

Previous theoretical studies of this phenomenon have 
used density functional theory ||^ which gives essentially 
the same results as Landau theory 0, namely a change 
from first to second order critical behavior with increas- 
ing potential strength. Continuous melting is allowed be- 
cause the external potential singles out density modes in 
one direction and the relevant Landau theory then con- 



tains only even powers in the local Fourier components of 
the density pg{^) = exp(iG-r). Unfortunately, the appli- 
cability of mean-field-like theories to problems with con- 
tinuous symmetry in 2d is limited since they drastically 
underestimate the effect of fiuctuations. Results from 
Monte-Carlo simulations are inconclusive: Although ear- 
lier simulations |10| claimed to have found a tricritical 



point at intermediate laser intensities and reentrance, re- 
cent studies from the same laboratory refute these 
results. We note that simulations with much larger num- 
bers of particles have still not completely resolved the 
nature of 2d melting without an external potential fl^ . 




FIG. 1. Triangular lattice with lattice constant a subject 
to a periodic potential (maxima indicated by dashed lines) 
for two different relative orientations: A) pAd — with 
o'a ~ V3a/2 and pa = 2, and B) psd = a'g with a'g = a/2 
and pb = 1. Also shown are low energy dislocations with 
Burgers vector b parallel to the corrugation of the potential. 

In this Letter we take a different approach, building 
on concepts developed in the context of dislocation me- 
diated melting theory |^,^. We model the experimental 
system by a continuum elastic free energy F — F^i + Fp, 
where F^i = ^ J (Pr[2^ufj + Au^^,] is the elastic energy 
associated with the colloidal displacement field u(a;, y) 
relative to the equilibrium position in the unconstrained 
solid characterized by the "bare" elastic constants /i and 
A. The quantity Uij — ^{diUj +djUi) is the 2d strain ma- 
trix. The effect of the laser potential with troughs run- 
ning along the x-axis is described by Fp = J (Pr U {y, Uy), 



where U{y,Uy) — —Uoa 



^{uy- 6- y)]. Uo mea- 



sures the strength of the laser potential and the mean 
colloidal spacing a is related to the particle density p by 
p — 2/-\/3a^. A similar model was used to discuss mod- 
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ulated superconducting films in Ref. Q . 

Let a! be the Bragg plane spacing for the orientation of 
the colloidal crystal (relative to the troughs) which pro- 
duces the lowest free energy (see Fig. |^). Incommensura- 
bility between a and the distance between the troughs, d, 
is accounted for by the mismatch parameter 6 = a' /d—p, 
where p is the integer closest to a' /d. 

For 5 = 0, we can write U{y, Uy) — —Uoa^^ cos[Kp • r], 
where Kp — ^pey is a. reciprocal lattice vector of the 
unperturbed system. If Uo = 0, there are algebraic Bragg 
peaks in the structure function of the crystalline phase 
and Mk = (Pk) ^ as the system size 

L^oo, where = '^J^^fxJ^^ l' I* straight- 
forward to show that A/k then vanishes for small Uq as 
A/k f/o Pk/(4-7;k) In contrast, A/k should al- 

ways vanish linearly with Uq in the liquid and hexatic 
phases [^,^ of the unperturbed colloid. The laser poten- 
tial will also induce long range bond orientational order 
in ipQ = (e^*^'''-') The bond order parameter ipQ van- 
ishes linearly with Uq in the liquid, vanishes like a power 
of Uq in the hexatic phase {ipe ~| Uq |S''6/(4-W6)^ where 
rjg is the exponent describing the algebraic decay of bond 
order), and approaches a nonzero constant as C/q — ^ in 
the solid phase. 

In our analysis of larger values of Uq , we shall also fo- 
cus on the commensurate case (6 — 0), where the spacing 
a' of the lattice planes parallel to the troughs equals an 
integer p times d. In Fig. |l], two particularly interesting 
orientations, denoted A and B, of the colloidal crystal and 
laser potential are shown. Only dislocations with Burg- 
ers parallel to the troughs have the usual logarithmically 
divergent energy. In orientation A, four of the six funda- 
mental Burgers vectors are disfavored by the potential, 
which requires that they be attached to a semi-infinite 
discommensuration string Q|. In orientation B, all six 
of the fundamental Burgers vectors are disfavored. The 
lowest energy Burgers vector parallel to the troughs has 
length -s/Sa. At sufficiently low temperatures, the laser 
potential is always relevant leaving massless phonon 
displacements Ux along the troughs with massive out-of- 
valley displacements Uy. We call this ordered phase a 
"locked floating solid" (LFS), reflecting its resistance to 
strains associated with Uy, and its ability to accomodate 
strains in u^- In reciprocal space, the LFS is character- 
ized by a structure function S{q) with a row of delta- 
function Bragg peaks at Gy""* = 2T:n/d [n e Z) along 
the q^-axis and power law Bragg peaks off this axis (see 
Fig. I). 

Upon integrating out the massive itj,-modes and using 
standard renormalization group methods [3|Jl5|] to elim- 
inate bound dislocation pairs in the LFS phase, we are 
left with a free energy with temperature and potential 
strength dependent effective elastic constants, 

^LFS = ^ / {K^s{dxUxf + ^i^a{dyUxf] . (1) 



Because only dislocations with Burgers vectors parallel 
to the troughs are important, the physics is described 
by an anisotropic scalar Coulomb gas or, equivalently, 
an anisotropic 2d XY model. For p smaller than a criti- 
cal value Pc we find that the LFS melts with increasing 
temperature at 



8^ 



(2) 



via an unbinding of dislocations with Burgers vectors 
b II Bj, before the laser potential which gives the mass to 
the Uj,-modes become irrelevant (sec below) . For T < 
the correlation function CG(r) = (Pg('')/^g(0)) recip- 
rocal lattice vectors with Gx shows power law decay 
CG{r) ~| (/icff/ifoff)i/2^2 + (^^g/^^g)i/2y2 |-^G/2_ The 

structure function is singular, S{q) ~ 1/ | q — G P^''*^, 
near the reciprocal lattice vector G The expo- 

nents {rjc,} are given in terms of the elastic constants 
by r/G = kBTGl/2Try^KcSlJ-oS- 
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FIG. 2. Schematic p = 1 phase diagram for orientations A 
and B. Insets: Schematic structure functions in the various 
phases. The x 's indicates delta-function Bragg peaks and 
circles algebraic peaks. 

Unlike conventional 2d melting |^, rye is universal at 
the melting transition, and given by ?7q = (G • b/47r)^, 
where b is the smallest allowed Burgers vector in the 
trough direction For orientation A, 6 = a and the 
exponent characterizing the algebraic order in the off- 
axis peaks (see insets of Fig. ^) closest to the gj,-axis is 
r^Q = 1/4; for the next row of peaks with Gx — An /a one 
gets ?7q = 1 , consistent with the algebraic decay observed 
in Ref. |^. For orientation B, 6 = \/3a and the six quasi 
Bragg peaks closest to the origin have different power 
laws; peaks on the x-axis have r/Q = 1, whereas the four 
off-axis peaks have rj^ — 1/4. Just above the melting 
temperature, we expect diverging translational correla- 
tion lengths (defined by the widths of Lorentzian peaks 
in the structure function) parallel and perpendicular to 
the troughs £,x ^ iy ^ exp[const./ | T — l^^"^]- 

To obtain the melting curve as a function of Uq we 
calculate K^g and ^^s and use Eq. (|^). We start from 
a microscopic model with a screened repulsive Coulomb 
interaction V{r) — Voo exp(— Kr)/r, where the screen- 
ing length is typically much smaller than a and Vq 
depends on the dielectric constant, k and the sphere ra- 
dius B. In such a dilute limit the two Lame coefficients 
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are equal, A = /i. Hence the Kosterlitz-Thouless melting 
temperature in the absence of a laser potential {Uq = 0) is 
simply given by ||^ fcsT'm = /ia^/67r. In the opposite limit 
of infinite potential strength the effective free energy sim- 
plifies to Eq. (|^) with Kcs = 3/i and ficS — A*- Using 



Eq. 



this gives keT^ = VSfia^/Sir « l.BfcijT". One 



might have thought that the melting temperature simply 
increases monotonically with Uq from to . It turns 
out, however, that fluctuations in the My-modes caused 
by lowering the potential strength from infinity typically 
increase the melting temperature. To see this, we in- 
tegrate out the Uj^-modes using the screened Coulomb 
potential; to leading order in /c^T/f/o a^nd Voe^^^/t/o we 
find (for orientation A) 
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23 + 104.A W;^ ^ 



where v = Voe-"'' /ksT, ^ = ^vkeTK^ and K = 3^*. 
Lowering the potential strength Uo always increases the 
shear modulus, whereas the behavior of the compres- 
sional modulus depends on the magnitude of v and na. 
When combined with Eq. (||), these expressions imply 
that the melting temperature Tm increases with decreas- 
ing Uo for KG > 5.6 (in Ref. § ko « 10), 
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thus implying reentrant melting for a band of temper- 
atures as a function of potential strength (see Fig. ||). 
For a more precise estimate of the phase boundary one 
would also need to determine the renormalization of the 
effective elastic constants by phonon nonlinearities and 
by bound dislocation pairs. In general, however, one 
expects only small downward renormalizations with in- 
creasing temperature IT^ ] which will not affect the exis- 
tence of reentrant melting. For small Uq, we find that the 
melting curve has a universal shape Tm{Uo) — Tm(0) 
[ln(fcBT„/[/o)]-i/^ with V w 0.36963 • • •. 

We now discuss the topology of the phase diagram for 
larger values of the integer commensurability parameter 
p. As illustrated in Fig. H, for p > pc a. floating solid 
(FS) with two soft phonon modes can intervene (bar- 
ring a direct first order transition) between the LFS and 
the liquid phase. To determine pc one must calculate the 
thermal renormalization of the laser potential, defined 
by Uois) = Uos^ exp[-^K^{ul)>] = Uos^", where \p 
is a renormalization group eigenvalue and the subscript 
">" indicates that the Uj^-modes are integrated over a 
momentum shell [A/s,A]. The laser potential becomes 
irrelevant at long wavelengths whenever Aj, < 0. The 
thermal average over the Uj,-modes requires an effective 
elastic long wavelength free energy appropriate to the 



FS phase. Because the laser potential breaks rotational 
symmetry, this coarse-grained free energy contains six in- 
dependent elastic constants, 

1 



+ 2-f0^ + 2a0uxy} , (6) 



:'^xxUxx ~r 2 yy^vy 



where 6 = ■^{dxUy — dyUx) is the local rotation angle 
induced by the phonon displacements. The resulting ex- 
pression for Ap is rather complicated, with an implicit 
dependence on the strength of the laser potential and the 
temperature However, a useful estimate results from 
neglecting the effect of the substrate on the elastic coef- 
ficients and considering an isotropic elastic free energy. 
For orientation A, one finds (| Ap = 2 - ^ p^ with 
the effective elastic constant fl — 2/i(2/i + X)/{3fj. + A). 
In the dilute limit ka 3> 1, relevant for many colloidal 
systems, we have fi = \ and the eigenvalue is negative 
for T > Tp, with j2 = 3/x/2 and keT^'^ « dua^/Airp^. 
When compared with the melting temperature of the LFS 
phase in the same limit of ko 3> 1 for weak potentials, 
ksT^ « /ia^/67r, wc find that a window of the FS phase 
exists for p > pc ~ 3^/3/2 « 3.7 for orientation A [ p^ . 
Note, however, that pc is not a universal constant but 
implicitly depends on the strength of the laser potential, 

Pc{Uo). 

For intermediate commensurate densities 1 < p < pc, 
no FS phase exists, and the dislocation unbinding tran- 
sition discussed above melts the LFS to a locked smectic 
(LSm) phase. In the LSm phase the laser potential is 
relevant and only one out of p possible troughs is pref- 
erentially occupied by the colloidal particles. At higher 
temperatures an Ising transition (p = 2) (more gener- 
ally, p-state clock model transition |1^]) takes the system 
to a modulated liquid phase. We expect the same kind 
of reentrant behavior for the LFS-LSm phase boundary 
as discussed above for p = 1. However, the LSm-Liquid 
phase boundary should depend more weakly on the po- 
tential strength: since the shear modulus is zero in the 
LSm phase the reentrance mechanism discussed above 
does not apply here. 

For p > Pc the laser potential becomes irrelevant before 
the LFS melts (via a roughening- like transition), leading 
to the uniaxially anisotropic FS described by the elastic 
free energy Fps, Eq. (||). Melting of this uniaxial solid 
(with 7 = a = 0) has been studied in detail by Ostlund 
and Halpcrin p9[ . Due to the uniaxial anisotropy there 
arc two different types of dislocations: type-I dislocations 
with Burgers vector along a reflection symmetry axis of 
the solid, and type-II dislocations at angles ±0o with 
respect to this axis. The laser potential in Fig. |i|A is ori- 
ented in such a way that type-I dislocations are parallel to 
the troughs of the laser potential, interact more weakly, 
and will therefore unbind first. This "type-I" melting 
transition leads to a "floating smectic" phase (FSm) sta- 
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bilized by the laser potential such that it retains quasi 
long range order in pc, = exp[iG • u] for G = 2'iTBylpd. 
We expect that a FSni phase intervenes between a FS 
and the liquid if p is sufficiently large {p > Pc = 4 |16[). 
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FIG. 3. Schematic phase diagrams for a) 1 < p < pc and 
h) p > p'^ > Pc (for orientation A). Heavy lines indicate 
phase transitions which are most likely of first order. Insets: 
Schematic structure functions. 

The novelty of this FS-FSm type-I melting is that, in 
contrast to the similar melting of the LFS, the destruc- 
tion of translational order in by dislocations takes 
place in the presence of a coupled "spectator" massless 
phonon mode Uy. To analyze this transition we proceed 
in the standard way Pp6| by introducing type-I disloca- 
tions into the elastic free energy and performing a duality 
transformation on the resulting Coulomb gas Hamilto- 
nian to convert it into a modified Sinc-Gordon model 



H 



\Vhl 
2K, 



- g cos [a{h + iilj)]\ + Ffs[4>\- (7) 



Here h(x^ y) is a dummy field which when integrated out 
in exp{—H/kBT) gives rise to a scalar Coulomb gas with 
interaction strength proportional to Kj. The quantity 
<p is the single valued part of the displacement fields, 
V'(q) = [Xxyql -ifi- j)ql)uyici)/q^ and g ^ 2e-^'l^-^ 
with the core energy of type-I dislocations. The melt- 
ing temperature can be determined from the condition 
that the renormalization group eigenvalue of g vanishes. 



S-t FS-FSm 



(8) 
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with suitable renormalized elastic constants. Here /i± = 
/i -|- 7 ± a and c is a dimensionless function of ratios of 
the FS elastic constants 15 1. The effect of the "specta- 



tor" phonon modes is to reduce the melting temperature. 
Just below Tps-FSm we find for the renormalized stiffness 

Kf{T) « iff (TFS-FSm)[l + const. (TpS-FSm " Tfl\ 

The coupling between Ux also induces cusps in the elas- 
tic constants of the spectator field Uy. Above the melt- 
ing temperature Ki is zero, type-I dislocations are un- 
bound, and at long length scales the elastic free en- 
ergy describing the floating smectic phase is of the form 
^FSm = \l (i^r{Bx{dxUyY + By{dyUyf} wlth Bx ~ 



Bx vanishes in the limit 7 = a = as one would expect 
for a rotationally invariant 2d smectic. The structure 
factor in the FSm phase exhibits power-law singularities 
at reciprocal lattice vectors G„ — 2Tmey/pd (for n ^ p), 
with an exponent jyFSm — ksT G1J2t:^ BxBy. The FSm 
phase melts via a second dislocation unbinding transition 
when Typsm > 1/4. 

We have only discussed laser potentials which are com- 
mensurate with the colloidal crystal. Sufficiently large 
incommensurability favors the discommensurations at- 
tached to Burgers vectors with an unfavorable orienta- 
tion relative to the troughs [Q , and leads to an enormous 
variety of phase transitions as a function of density at 
fixed trough spacing which we hope to discuss in a fu- 
ture publication. 
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